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ABSTRACT 

In this paper we examine the profiles of the complex Na i D-lines in the Red Rect- 
angle. The spectra were acquired with the ARCES echelle spectrograph {R = 38, 000) 
on the 3.5-m telescope at the Apache Point Observatory. Additional spectra taken with 
STIS were acquired from the Hubble Legacy Archive (HLA) and were used to indepen- 
dently confirm the spatial origin of the spectral features. The profile of a single D-line 
consists of double-peaked emission, red-shifted absorption and blue-shifted absorption. 
We find that the double-peaked emission originates from the bipolar outflow, the red- 
shifted absorption feature is due to the photospheric line, and the blue-shifted absorp- 
tion arises from the bipolar outflow as seen against the photosphere of the luminous 
post-AGB component in HD 44179. In order to better understand the Nai D-line pro- 
file, we examined the periodically variable asymmetric photospheric absorption lines. 
The asymmetric lines are interpreted as a signature of slow self-accretion following 
enhanced mass-loss around periastron. An empirical model was constructed to remove 
the photospheric component from the Na i D-line profile in order to study the nebular 
emission and absorption of sodium along the line-of-sight to the primary. This paper 
also discusses the different origins of the single-peaked emission, the double-peaked 
emission and the blue-shifted and red-shifted absorption components. 

Key words: ISM: jets and outflows ~ stars: AGB and post-AGB - stars: individ- 
ual(HD 44179) - stars: mass-loss - binaries: close 



INTRODUCTION We will refer to th is optically-thick disc, resolved by HST 

( Co hen et al.ll2004n and grou nd-based interferometric obser- 



The Red Rectangle (RR) , an enigmatic protoplanetary neb- yations jRoddier JID M : lOsterbart. Langer. fc Weieeld 



ula, reveals itself to be a complex astrophysical system. The |,„„4, ,, . ,. i- m, i- ■ i 

' ^ , , , ,^^ ..,„„.,. . I19971 I. as tiie circumbmary disc, ine dimensions assumed 

central source oi the nebula, HD 44179, IS a binary system m r ^^ ■ i- i- ^ i i- ^ r-nn 

■' ■' tor the circumbmary disc are based on a distance ot 710 pc 

a short-lived evolutionary phase between the asymptotic gi- 177 — n — T^ ; — rllr>,^/^rJ\ Tin -i .^ ^■ , 

, /.^T^N , , , , , ^. ,.t (IMen'shchikov et al.ll2002l ). While the distance i s uncertain, 

ant branch (AGB) and the planetary nebula stage 01 its lite, 1 . .1 • 1 . • ,. ,■ n • — \ — '1 — \ — 71 

, / , , , , , we adopt this value m our investigation. Buiarrabal et al.l 

when the expelled circumstellar matter IS ionized by the pho- ,r,nr,tr\ i j^t, ^ j- r^i • u- j-i 

. , , I — ,,,. , , „■' „ „„„„ (2005) traced the outer radius ot the circumbmary disc to 

tometric primary stars hot core (Ivan Winckcl 2003, 2007: , „_„ .tt rm j- r ^i ^ 1 -^ ■ ^i • 1 ■ 

U. , , — TT"^ 7\\ ~„A r^ 7^ "rn U In r„, , . 1,850 AU. 1 he radius 01 the central cavity m the circumbi- 

ISiodmiak et al.ll2008l : ISzczerba et al.l 120071 ). The dominant j- ■ u^ • j i- 4-4- j r 1,^ j w u 

~ ~~ — ^~ — T-, r \ , — , , . . , nary di sc is obtained trom a sc attered light model to be 

teatures ot the nebula s torm are the X -shaped bicomcal out- .. .TiiTT — , , n.., : — ninnnoK mu it.- 1 r -l-l. ■ 

, , „, , , r-— ; n uJT^A N , 14 AU HMen shchikov et al.ll2002l l. The thickness of the cir- 

nows nearly 2 on the sky (lUohen et al.ll2004l l I that emerge ,. .- , , , ^ • 1 ^ 1 1 ^ r.^ att 

^ ^ " cumbmary disc has been d etermined to b e between 90 AU 



from a disc (roughly 0. 1" on the sky) , seen nearly edge-on. JMen'shchikov et all B a nd 107 AU jBuiarrabal et al 

20051: iRoddier et all 1 19951: lOsterbart. Langer. fc Weigeh 
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1997]). It is worth noting that values of the viewing angle 
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vary from ~ 0° to 11° JLopez. Mekariiia. fc Lefevrelll995l : 
iRoddier et"aIlll995l : lMen'shchikov et al.ll2002h . These view- 
ing angles would correspond t o incli nation angles between ^ 
90° to 79°. lBuiarrabal et alj (|2005l ') find that an intermedi- 
ate value of 5° fits their data best, which would correspond 
to an inclination angle of 85° . 

At the centre of the cavity in the circumbinary disc 
are the components of the binary system: the photomet- 
ric primary (a mass-losing post-AGB star), and the sec- 
ondary (a main sequence star with an accretion disc and 
high- velocity jet). For simplicity we shall refer to the pho- 
tometric primary as the primary star throughout this pa- 
per. The primary star has an effective temperature of 
7,700 K, an est imated mass of 0.8 M^, a luminosity of 
about 6.0 00 Lm IIMen'shchikov et al.ll2002l l. and a radius of 
0.21 AU (I Witt et al.ll2009l ). The secondary has a mass of 
roughly 0.9 Mq a nd is believed t o be a main seque nce star 
jWitt et aLlbOOgJ ). In the work of IWitt et all (|20()9l 'l the de- 
pendence of the Ha line on orbital phase was investigated, 
revealing a high velocity jet {vmax ^ 560 km s~^) originat- 
ing from an accretion disc around the secondary. The accre- 
tion disc around the secondary is powered by mass transfer 
from the primary sta r via Roche-lobe overflow at perias- 
tron (|Witt et al.ll2009l 'l. The ionizing photons from the hot 
inner regions of the accretion disc provide the necessary en- 
ergy to excite the Extended Red Emission (ERE). The ERE 
requires far-ultra violet photons with energies E > 10.5 eV 
l|Witt et al.l 120061 ). The accretion disc also provides the en- 
ergy source to maintain the Hll region inside the cavity of 
the optically-thick circumbin ary disc that was suggested by 
IJura. Turner, fc Balml l|l997l ). 

The near edge-on optically-thick circumbinary disc pre- 
vents a direct line-of-sight view of the binary. Light origi- 
nating inside the cavity can only be observed via scattering 
off material located above and below the disc. This obscured 
view of the binary, illustrated in Fig. [T] will be referred to 
as the indirect line-of-sight and is responsible for the difB- 
cult interpretation of the spectra. The indir ect line-of-sight 
to the central object was first studied by IWaelkens et al.l 
l|l996h , who determined the effective inclination angle to be 
35° . The effective inclination angle is the angle at which light 
from the central binary is scattered into our line-of-sight by 
the material near the top and bottom of the disc. This view- 
ing geometry has been confirmed by high-resolution op tical 



l|Osterbart. Langer. fc Weigelt 



infrared imaging ( Tuth ill et al 



sd by high-resolution op tical 
1997 ^ : ICohen et al.ll2004l ) and 



20021 ). For light from the cen- 



tral source (which can only been seen via scattering), we 
adopt the effective inclination angle of 35° in our calcula- 
tions. 

The outer regions of the nebula are, however, seen di- 
rectly. To properly study emission processes that originate 
in the outer nebula one must use the inclination angle of 
the nebula (85°) rather than the effective inclination angle 
(35°). 

Since HD 44179 can only be studied via scattered light, 
it is important to discuss the known uncertainties. The ef- 
fective angle at which we are viewing the stars is likely the 
mean scattering angle. The variation of the scattering angle 
over the duration of the orbit, which best fits the 0.14 ma g 
variation in the light curve, is 1.6° ( Waelkens et al.lll996l ). 
The average value for the scattering angle corresponds to 
the effective inclination angle of 35° . In addition, the tilt of 



the system implies that the scattering angle for light from 
the north and the south of the nebula is slightly different 
(of order 5°). Therefore, the effective inclination angle is 
just that, an effective value that we use in order to make 
progress in understanding this complex system. The varia- 
tion in the value for the effective inclination angle will cause 
uncertainties in; the orbital elements, the derived velocities, 
and the mass of the secondary. 

The properties of the stars also depend upon t he es- 
timation of the luminosity (|Men'shchikov et al.ll2002l ). The 
luminosity depends on the uncertain distance to the star, 
assumptions about the geometry of the circumbinary-disc, 
and the evolutionary state of the system. Since this a close 
binary, it is likely that there has been some interaction in the 
past resulting in a complex evolutionary history. Therefore, 
it is likely that the evolutiona ry track of the star s differ 
fro m that of single stars, as in IVassiliadis fc Wood! (|l994r ) 
and [Bloeckeg ([1995|). As a result, the quoted values that we 
have used for the mass and luminosity of the primary re- 
main uncertain. At present, these uncertainties can not be 
addressed with current observational data. 

Studying the nebular outflow and its mechanisms is im- 
portant for understanding how stellar material is returned 
to the interstellar medium. The detectable emission of the 
Ha recombination line in our spectra is limited to the inner 
parts of the nebula. The Ha emission is not observed in the 
bi-conical lobes far away from the central source. Studying 
the Nal D-lines, easily excited resonance lines, allows us to 
probe the mass-loss to much greater distances from the cen- 
tral source. The Nal D-lines also allow us to directly study 
the mass-loss from the primary star. 

The details of the observations, and the radial velocity 
(RV) measurements are given in § [2] The determination of 
the orbital parameters is discussed in § O In § U we discuss 
the Nal abundance and the model atmosphere used to de- 
termine the photospheric line profiles. The behaviour of the 
photospheric lines and the processing of light in the RR is 
discussed in § [5] The origins and measurements of the com- 
ponents of the Na I D-line profile are discussed in § |6l In § [7] 
the origin of the single and double peaks are discussed. The 
conclusions are presented in §|8] 



2 OBSERVATIONS AND MEASUREMENT 
2.1 Ground-based observations 

2.1.1 Apache Point Observatory 

In the present paper we examine 33 high-resolution echelle 
spectra of HD 441 79. The data set includes the 17 spectra 
used in the work of iWitt et al.l l|2009r) . The observation dates 
are detailed in Table [T] The orbital coverage of the obser- 
vations is illustrated in Fig. (2] The details of Fig. [2] will be 
discussed in §[3] The spectra were obtained with the ARCES 
echelle spectrograph ([Wang ot al. 2003) on the 3.5-m tele- 
scope at the Apache Point Observatory. The spectral range 
of the spectrograph is 3,700 A to 10,000 A at a resolving 
power oi R — 38,000. The velocity resolution is 8 km s^^. 
The entrance slit used measures 1.6" x 3.2" on the sky and 
cannot be rotated. The relative size of the slit compared to 
the nebula can be seen in Fig. [S] The first spectrum in the 
combined data set was acquired on 2001 February 6. The 
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Figure 1. This to-scale diagram of the inner region of the RR 
was constructed for the distance of 710 pc. The outer radius of 
the circumbinary disc, 1,850 AU, is not shown. The central binary 
is displayed a factor of 10 larger than the scale for the rest of the 
system. The inclination angle, angle between the plane of the sky 
and the equatorial plane, is 85° . The direct and indirect lines-of- 
sight are also indicated for both the northern and southern ends 
of the nebula. The cavity is filled with the low- velocity, mass-loss 
driven outflow of the primary star and the high-velocity jet from 
the accretion disc of the secondary. The location of the primary 
star is shown for the orbital phase of inferior conjunction; see 
Fig. [2] for a to-scale diagram of the binary orbits depicting the 
orbital phases. 



most recent spectrum was acquired on 2009 April 7. The 
firs t 7 and 17 sp ectra, respecti vely, were studied previously 
by iHobbs et af] (|2004i ) and Witt et al.l 12009 ). Almost all 
the new spectra were taken by DGY, based on prescribed 
d ates from ANW, r eferring to the results and implications 
of I Witt et ai] l|2009l ). The data in Table [l] are based on re- 
duction of the raw spectra to normalized spectra by JAD 
and measurements for individual nights by LMH. The reduc- 
tion of the spectra and inst rumental set-up are un changed 
since the previous work by iThorburn et al.l (|2003l ) and by 
iHobbs et ai] (|2004l '). 

The present study focuses on the complex Nal D-lines, 
as exemplified in Fig. [l] The large number of observations 
were needed in order to study variability in the spectra over 
the course of the binary orbit, not for co-adding. Each ob- 
servation is comprised of 2 averaged spectra, with signal-to- 
noise (S/N) at the continuum in the region of the Nal D-lines 
of order 300 for the combined spectrum. We have also exam- 
ined several absorption lines with no emission component, 
Fig- El In our ground-based spectra, the slit (1.6" x3. 2") re- 
ceives both the scattered light from the primary star and 
the emission from the inner nebula sinmltaneously. 



2.1.2 McDonald Observatory 

A spectrum of the Nal D-lines was also acquired with the 
Harlan J. Smith 2.7-m telescope at the McDonald Obser va- 
tory (McD) using the TuU spectrograph (|Tu11 et al.lll995l ') in 
its high-resolution mode (ts21). The observation was made 
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Figure 2. Based on the orbital parameters shown in the column 
M2 of Table [2] the orbits of the binary components have been 
plotted to-scale (assuming a distance of 710 pc) and are shown 
from the privileged vantage point above the orbit (inclination i = 
0) . The solid ellipse indicates the orbit of the primary, while the 
dashed ellipse indicates the orbit of the secondary. The hash at the 
origin marks the centre-of-mass. Orbital phases are designated on 
the primary orbit with the major and minor hashes. The direction 
of motion for the primary star is counter clockwise. Apastron (A) 
<j> = 0.79, and periastron (P) {j> = 0.29 occur at the intersection 
of the orbit and the line passing through the antipodes of both 
orbits. The direction to Earth is to the right and 5° below the 
plane of the page. Inferior conjunction (IC) (p = 0.21, and superior 
conjunction (SC) </< = 0.59 occur at the intersection of the arrow 
indicating the direction to Earth and the orbit. The terminology 
of IC and SC was chosen to be consistent with IWaelkens et al.l 
(1996) and lWitt et al.l l|200^ . Maximum red-shift of the primary 
occu rs at (f> = 0.42. The loc ation of phase zero was chosen to agree 
with lWaelkens et al.l {idM) and corresponds to JD 2448300. The 
circles on the orbits indicate the locations of the stars for all 33 
APO observations. All distances are in astronomical units (AU). 
Note that the diameter of the post-AGB primary is estimated to 
be 0.42 AU. 



on 2010 December 24. The crossed-dispersed echelle spec- 
trometer was configured with the 79 g mm~^ grating (El), 
the 145 ^m slit (Slit 2), and a 2048 x 2048 CCD (TK3). 
The slit is approximately 0.34" x 8. 2" on the sky. The slit 
orientation on the sky is not fixed. Order 39 of the El grat- 
ing spectrum was centred at 5890 A. The resolving power 
(R « 200,000) was estimated by measuring the FWHM of 
the unblended lines in the Th-Ar comparison spectrum. The 
corresponding velocity resolution is 1.5 km s~^. The orbital 
phase of the observation was = 0.814. The data were re- 
duced and wavelength-calibrated with standard iraf rou- 
tines. Plotted in Fig. [6] is the averaged profile for the Nal 
5889.950 A line, with a continuum S/N ratio of approxi- 
mately 50 near the sodium lines. The average profile was 
constructed from four 1800 s exposures. 
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Figure 4. Tlie spectra siiown were taken at APO witii a resolving power of ij = 38,000. Panel a shows the Nal D2 5889.950 A line, 
while panel b shows the Dl 5895.924 A line. The spectra are the observed, averaged, and unsmoothed data with S/N of roughly 300. 
The spectra are shown at the phases of </> = 0.423 (solid line) and = 0.904 (dashed line) corresponding to the maximum and minimum 
red-shift, respectively, of the absorption associated with the photosphere of the primary star. See Fig. [2] for the respective locations of 
the stars in the binary system during these phases. The flux is in arbitrary units (a.u.). The heliocentric velocity can be determined by 
adding 19 km s~^ to the centre-of-mass velocity. 
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Figure 6. The profiles shown were acquired at McDonald Observatory with a resolving power R ^ 200,000. Panel a shows the Nal D2 
5889.950 A line, while panel b shows the Dl 5895.924 A line. The spectra are the observed, averaged, unsmoothed data and have a S/N 
of approximately 50. The spectra are shown at the phase of </> = 0.814 (Fig. |2]l. The flux is in arbitrary units (a.u.) . The heliocentric 
velocity can be determined by adding 19 km s~^ to the centre-of-mass velocity. 



2.2 Space-based observations 

2.2.1 HST/STIS 

Nal D spectra obtained with STIS can be seen in Fig. [T] 
These spectra were obtained through the Hubble Legacy 
Archive. The spectral resolution of the STIS spectra is lower 



than that of the other spectra presented in this paper, but 
they have the superb spatial resolution of the HST along the 
slit. The 1440 s exposure of the Nal D-lines was taken on 
1999 February 12 by Hans van Winckel (proposal ID 7297). 
The 52"x0.2" Fl aperture and the G750M grating with a 
resolving power oi R = 5, 000 were used to acquire this spec- 
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Table 1. Details of the APO observations and measurements 



# 


Year 


Month 


Day 


UT 


HJD 


AHJD 


Phased 


AyheUct 


AVlsa 


RV^ 


s^y stronc^ 


RV^ea*;^ 


1 


2001 


2 


6 


4:21:46 


2451946.68531 





1.470 


-16.4 


-35.3 


27.2 


30.6 


13.7 


2 


2002 


12 


31 


7:15:39 


2452639.80723 


693.12 


3.650 


-1.6 


-20.4 


17.3 


14.4 


23.2 


3 


2003 


11 


24 


8:25:03 


2452967.85463 


1021.17 


4.681 


13.9 


-4.9 


17.3 


14.4 


23.1 


4 


2004 


1 


3 


8:15:46 


2453007.84896 


1061.16 


4.807 


-3.0 


-21.8 


11.9 


9.4 


16.7 


5 


2004 


1 


10 


4:28:03 


2453014.69072 


1068.01 


4.829 


-5.6 


-24.4 


12.3 


8.3 


20.9 


6 


2004 


1 


31 


4:58:31 


2453035.71118 


1089.03 


4.895 


-14.2 


-33.0 


10.7 


8.1 


16.1 


7 


2004 


2 


3 


4:21:35 


2453038.68538 


1092 


4.904 


-15.2 


-34.0 


10.7 


6.5 


23.3 


8 


2006 


11 


12 


10:33:59 


2454051.94356 


2105.26 


8.090 


17.7 


-1.2 


13.9 


12.4 


21.1 


9 


2006 


12 


28 


5:36:38 


2454097.73848 


2151.05 


8.234 


-0.1 


-18.9 


20.4 


23.0 


9.5 


10 


2007 


1 


15 


5:49:11 


2454115.74694 


2169.06 


8.291 


-8.0 


-26.8 


25.3 


28.3 


10.3 


11 


2007 


3 


28 


1:58:16 


2454187.58204 


2240.9 


8.517 


-24.8 


-43.7 


26.0 


31.3 


12.0 


12 


2007 


10 


3 


11:45:44 


2454376.99043 


2430.31 


9.112 


24.7 


5.8 


14.2 


13.3 


20.6 


13 


2007 


11 


5 


8:42:34 


2454409.86573 


2463.18 


9.216 


19.9 


1.1 


18.1 


17.8 


17.8 


14 


2007 


12 


18 


8:22:37 


2454452.85368 


2506.17 


9.351 


4.1 


-14.7 


29.6 


33.9 


11.1 


15 


2008 


1 


10 


5:30:43 


2454475.73424 


2529.05 


9.423 


-5.7 


-24.6 


31.0 


36.4 


12.9 


16 


2008 


2 


8 


6:51:57 


2454504.78953 


2558.1 


9.514 


-17.1 


-36.0 


26.2 


31.4 


13.1 


17 


2008 


2 


29 


2:18:08 


2454525.59801 


2578.91 


9.580 


-22.2 


-41.0 


21.3 


18.7 


27.1 


18 


2008 


12 


16 


7:15:49 


2454816.80727 


2870.12 


10.496 


4.8 


-14.0 


26.0 


31.8 


10.5 


19 


2008 


12 


19 


7:34:18 


2454819.82015 


2873.13 


10.505 


3.4 


-15.4 


24.6 


31.8 


10.5 


20 


2008 


12 


25 


5:59:41 


2454825.75448 


2879.07 


10.524 


1.0 


-17.9 


24.6 


30.7 


8.0 


21 


2008 


12 


28 


5:48:06 


2454828.74644 


2882.06 


10.533 


-0.4 


-19.2 


24.8 


29.5 


13.1 


22 


2008 


12 


31 


6:39:41 


2454831.78225 


2885.1 


10.543 


-1.8 


-20.6 


23.3 


28.6 


12.0 


23 


2009 


1 


2 


6:28:38 


2454833.77109 


2887.09 


10.549 


-2.7 


-21.5 


22.8 


11.2 


18.6 


24 


2009 


1 


5 


7:04:07 


2454836.79917 


2890.11 


10.558 


-4.1 


-22.9 


23.1 


20.7 


28.9 


25 


2009 


1 


7 


6:26:23 


2454838.77271 


2892.09 


10.565 


-4.9 


-23.7 


23.2 


20.7 


29.3 


26 


2009 


1 


10 


6:17:57 


2454841.76703 


2895.08 


10.574 


-6.2 


-25.0 


22.9 


21.9 


28.9 


27 


2009 


1 


14 


6:47:14 


2454845.78727 


2899.1 


10.587 


-7.9 


-26.8 


20.8 


19.7 


26.4 


28 


2009 


1 


17 


6:33:11 


2454848.77743 


2902.09 


10.596 


-9.2 


-28.0 


19.8 


18.2 


24.4 


29 


2009 


1 


30 


2:14:41 


2454861.59742 


2914.91 


10.636 


-13.8 


-32.6 


18.3 


16.8 


30.5 


30 


2009 


2 


2 


1:51:48 


2454864.58138 


2917.9 


10.646 


-14.8 


-33.7 


18.2 


15.6 


23.4 


31 


2009 


2 


6 


1:25:21 


2454868.56281 


2921.88 


10.658 


-16.2 


-35.0 


17.0 


15.6 


22.8 


32 


2009 


2 


13 


2:40:55 


2454875.61487 


2928.93 


10.680 


-18.5 


-37.3 


17.6 


15.0 


21.9 


33 


2009 


4 


7 


2:00:20 


2454928.58261 


2981.9 


10.847 


-24.4 


-43.2 


12.7 


7.3 


25.3 



"^ Based on a period of 318 days; the tabulated phase values are orbit number before the decimal point, followed by the 
phase value. The orbit numbers start with this observing campaign. Phase corresponds to JD 2448300. Refer to Fig. [2] 

" All velocities are in km s~^. 

'^ RV determined as in lHobbs et al.l l l2004h . 

RV measured for the strong component of the C I 5052 A line. 
^ RV measured for the weak component of the C I 5052 A line. 



trum. The velocity resolution of these data is 60 km s~^. 
The S/N at the continuum for both the sodium and hy- 
drogen spectrum presented in the present paper is of order 
200. The occuhing bar measures 0.5" x 0.2". The length 
of the slit is 52". This coronagraphic slit spectrum was ori- 
ented along the north-east whisker, position angle 42.056°. 
This slit orientation is similar to one of the slit orienta- 
tions in ISchmidt fcwitti l|l99lh . Fig. [3] is a to-scale dia- 
gram illustrating the orientation and location of the vari- 
ous s lits as projected on to the nebula. The relevant slits 
from ISchmidt fc WittI (jl99ll ). the STIS observations, and 
the relative size of the APO slit and of the McD slit are all 
illustrated in the Fig. |3] The nature of the coronagraphic 
spectra allows one to study the spatial origin of the Na I D- 
line components. The STIS spectrum was acquired at the 
orbital phase tji = 0.19. 

The 2029 s exposure of Ha, shown in Fig. |8] was taken 
on 1998 March 26 by Theodore GuU (proposal ID 7593). 



The spectrum was acquired with the same aperture, grat- 
ing and position angle as the Nal spectrum. This spectrum 
was also obtained from the Hubble Legacy Archive. This 
observation was made during the orbital phase tj) = 0.17. 
In all STIS spectra the northern and southern regions of 
the slit spectrum were extracted as separate apertures with 
standard iraf routines. The region very close to the coron- 
agraphic bar was avoided in the aperture selection, due to 
contamination of the spectrum from scattering by the coro- 
nagraphic bar. The aperture selection for north and south 
apertures started at about l" away from the centre of the 
occulting bar. 



2.3 Radial velocity 

The stellar radial velocity (RV) for each observation is re- 
quired in order to remove the photospheric absorption line 
from the observed profile. The removal of the photospheric 
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Figure 3. A to-scale diagram showing the RR as seen on the sky. 
The STIS slit is to-scale length wise, but not width wise, being 
only 0.2" wide. The width of the coronagraphic bar i s to-scale. 
The two MMT slit locations of ISchmidt fc WittI l|l99ll '). that are 
mentioned in this paper, are shown to-scale. The size of the APO 
slit is shown on the lower right. The McDonald (McD) slit, also 
shown on the lower right, is to-scale length wise, but not width 
wise owing to its small width on the sky (0.34"). Both the APO 
and McD slits were centred on HD 44179 during observations, the 
orientation of these slits on the sky cannot be controlled. 
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Figure 7. This spectrum of the Nal D-lincs was taken with STIS 
and has a resolving power ij = 5, 000. The north and south aper- 
tures were extracted separately. The centre-of-mass velocity scale 
is centred on the Nal D2 5889.950 A line. To convert back to he- 
liocentric velocity, add 19 km s~^ . These spectra were acquired 
during the phase = 0.19, see Fig. [21 
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Figure 5. Shown is a comparison of three observed, averaged, 
unsmoothed, photospheric lines. These data were acquired with 
APO and have a resolving power of R = 38, 000. The C I 
5052.167 A profiles were divided by 1.5 in order to show them 
on the same vertical scale. Note the similarity between different 
species, ionization states and wavelengths. On the far right, each 
profile is labelled with observation number, while on the far left 
the phase is shown. Spectra # 13 and # 28 occur near inferior 
and superior conjunction respectively. Spectrum # 7 occurs near 
minimum red-shift of the primary and # 15 occurs at maximum 
red-shift. Fig. [2] illustrates the locations of these key orbital po- 
sitions. To convert back to heliocentric velocity, add 19 km s^^ 
(the centre-of-mass velocity of the system) . 



Figure 8. Shown are the STIS Ha line profiles for the north and 
south apertures extracted separately. The resolving power of the 
spectra is i? = 5, 000. To convert back to heliocentric velocity, 
add 19 km s~^ . These spectra were taken at phase 4> = 0.17, see 
Fig.H 



absorption component from the Nal D-line profile is dis- 
cussed in § 15.51 The RV also allows us to determine the or- 
bital parameters of the system. The stellar RV of each obser- 
vation was measured in two ways. Method one (Ml), as de- 
scribed below, involves measuring the centroids of three se- 
lected lines in the spectrum of the primary. Ml is less prone 
to peculiarities of a single line as those would be averaged 
out. However, Ml is affected by the periodically asymmet- 
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ric wing present on these absorption lines. Method two (M2) 
accounts for the line asymmetries by using a two-component 
Gaussian to deconvolve one of these lines. Possible causes of 
the asymmetry are discussed in detail in §(5] 

Th e first method (Ml ) was identical to the method out- 
lined in I Witt et al.l (|2009l ). In order to construct a complete 
and consistent dataset for all 33 observed spectra in the APO 
dataset, we chose to continue using Ml. The centroid of each 
line, measured with iraf, was used to determine the RV. 
The three absorption lines used in those measurements were: 
Fell 4923.927 A, Ci 4932.049 A, and Ci 5052.167 A, and 
the RV results hav e typical random errors of ± 0.4 km s~^ 
iJHobbs et al.ll2004l ). These lines were chosen because they 
are not too weak to be measured, nor are they too near sat- 
uration. These lines are also not contaminated by telluric 
features, they are not blends with lines from other elements, 
and are not unresolved blends with other multiplet lines. 
The RV reported is the average of the RV determined for 
each of these lines. These values are tabulated in Table [T] 
under the column heading RV, and are plotted in Fig. [9] 

The profile asymmetry affects the determination of the 
RV when using the entire profile, as in Ml. The second mea- 
surement technique (M2) uses a two-Gaussian deblending 
of the Ci 5052.167 A line to separate the strong line core 
and the weaker wing. A sample of the line profiles, shown 
in Fig. O reveals that the asymmetric profiles of nearly 
all photospheric absorption lines vary as a function of or- 
bital phase in the same manner. We recognize that a two- 
component deblending is simplistic; however, it does allow 
us to quantify the asymmetric lines. This C I line was cho- 
sen as representative of the sample of lines and met the 
criteria listed in the previous paragraph. Since the lines are 
asymmetric during most of the orbit, this has the effect of 
damping the RV amplitude at maximum and minimum RV. 
This exp lains the difference between Ml RV and RV mea- 
sured bv lvan Winckel. Waelkens. fc Waters! l| 19951 ) at maxi- 
mum and minimum RV, as shown in Fig. |9l The FWHM of 
the Gaussian for both profile components was assumed to be 
the same, assuming the primary star is the dominant source 
of light in the central cavity. The width of the symmetric 
profile spectrum # 13, which occurs at inferior conjunction 
{(f) = 0.21), was used as the intrinsic width of the underly- 
ing line profile. From our observations, we measure a nearly 
constant total equivalent width having no dependence on the 
orbital motion of the star. Therefore, we also constrained the 
equivalent width of the line to be constant. We associate the 
strong component with the true RV of the primary star. 

The resulting strong (RVstrong) and weak (RV^eafc) 
components are tabulated in Table [l] The RVstrong val- 
ues are plotted in Fig. [9] (filled triangles). The max- 
imum R V values of M2 are in agreem ent with the 

data of Ivan Winckel. Waelkens. fc Waters! IJ1995! ). while 
the minimum RV values of Ml ( grey circles) match the 
Ivan Winckel. Waelkens. fc Water3 (|l995l ) data better. This 
indicates that M2 is at least reasonable for determining the 
RV of the primary star. The results of M2 for the strong 
component are equivalent to measuring the RV at the deep- 
est part of the absorption line. The velocity of the RVmeafc 
with respect to the RVstrong are plotted in panel c of Fig. [10] 
and will be discussed in § 15.51 




1.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
Phase 



Figure 9. Radial velocity curves for HD 44179, comparing the 
differ ent measurements techniques. The open squares are data 
from Ivan Winckel, Waelkens. &: Waters! lll995l ), the grey circles 
are m easurements made in the same manner as iHobbs et al.l 
J2004|) , using Ml. The filled triangles represent the deblended 
strong Ci component measured with M2. See Fig.|2]for the loca- 
tion of orbital phases. 



3 ORBITAL ANALYSIS 



A Lomb-Scargle periodogram ijScargld Il982h was con- 
structed for all 33 BVstrong measurements presented in this 
paper. The resulting period found was 317 ± 2 days, which 
is in agreement with the 318 ± 3 day period measured by 
IWaelkens et al.l (|l996! ) . This period was used as a c onstraint 
in the VCURVE program (JBertiau fc Grobber i 1969). 

The orbital parameters for the binary system were de- 
termined with the RV values for both Ml and M2. Ml uti- 
lized the RV measurements listed in Table [T] and were fit 
using a modified version of the FORTRAN program vcuRVE. 
M2 used VCURVE to fit the strong component (RVstrong). 
The orbital parameters determined from vcuRVE for both 
methods are shown in Table [2l The ve locity values from 
Ivan Winckel. Waelkens. fc Waters! (|l995l ) are also presented 
in Table [5] The quantities below the horizontal rule (ai, 
Mt, and m2) in Table [2] are derived using the effective in- 
clination angle of 35° and a primary mass of 0.8 Mq. Aside 
from the orbital period difference, the orbital parameters 
found via M2 are in very cl ose agreement with those of 
Ivan Winckel. Waelkens. fc Wa ters (1995). This close agree- 
ment suggests that the RV measurements using M2 pro- 
duce parameter s for the system are more consistent with 
IWaelkens et al.l (|1996! ). It should be noted that the uncer- 
tainty in the mass function and secondary mass are most 
sensitive to the uncertainty in the semi-amplitude. 

Based on the orbital data determined using the vcuRVE 
fit to the strong component (RVstrong), a to-scale schematic 
of the binary orbits has been constructed to aid in the in- 
terpretation of the spectra, see Fig. (2] The locations of the 
33 APO observations to date are also indicated in the figure 
along with key orbital locations. Orbital phase zero, = 0, 
was chosen to correspond with the location of phase zero 
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Figure 10. In panel a, the Roche lobe radius in astronomical 
units (Rij) is plotted as a function of orbital phase as calcu - 
lated by equation (46) of lSepinskv. Willems. fc Kalogeral ll2007f) . 
In panel b, the derivative (dRij dt~^) of the curve in panel a is 
shown. Positive slope indicates an expanding Roche lobe radius, 
while a negative slope indicates a contracting Roche lobe radius. 
In panel c, the velocity of the weak component with respect to 
the primary star (AV) is plotted so that the negative velocities 
are above the zero-line. Negative velocities indicate material mov- 
ing away from the primary star and vice versa. The vertical lines 
represent the phases of inferior conjunction (IC) </< = 0.21, peri- 
astron (P) (fi = 0.29, maximum RV (maxRV) </< = 0.42, superior 
conjunction (SC) 4> = 0.59 and apastron (A) <f> = 0.79. See Fig. [2] 
for the location of orbital phases. 



defin ed by IWaelkens et al.l (119961 ) to aid in comparison of 
data. IWaellcens et al.l ( 19961 ) chose 4> — to correspond to 
JD 2448300, near minimum RV. 

Since the components of the binary are in an eccen- 
tric orbit, the Roche lobe radius is not constant as a func- 
tion of orbital phase. The consequences of the variable size 
of the Roche lobe are discussed in § 15.31 In panel a of 
Fig. \TU\ the Roche lobe radius of the primary star was cal- 
culated as a function of orbital p hase using equation (46) of 
ISepi nskv. Willems. fc Kalogeral (|2007l ). The exact value of 
the effective Roche radius may differ, as will be discussed 
later; however, the relative variation and orbital phase de- 
pendence should remain the same. 



PHOTOSPHERIC ABUNDANCE OF 
SODIUM 



In order to remove the photospheric absorption line from 
the Nal D-line profile (see Fig. |4} the Na abundance is re- 
quired for the construction of a photospheric model line. 
The complex line profile of the Nal D-lines were not suit- 
able for the determination of the Na abundance. The near in- 
frared (IR) photospheric absorption multiplet at 8183.256 A 
and 8194.824 A provided an alternate approach in the de- 
termination of the Na-abundance. The spectral coverage of 
ARCES included the 8183.256 A and 8194.824 A Nal lines; 
however, the 8183.256 A line was buried in a forest of telluric 
features. Fortunately, the 8194.824 A line, the stronger of the 
two, fell between the telluric features for some of the obser- 
vations and we were able to measure its equivalent width to 
be 80.7 ± 3.0 mA. 



The PHOENIX code, Version 15.04.00E llHauschildt 
1992| . Il993l: [Hauschildt fc Baronlll995l: lAllard fc Hauschildt 



I995I: iBaron et al.l 



Hauschildt. Baron, fc Allard 



199a iHauschildt et al.l Il996l: 
I997I: iBaron fc Hauschildt 



1998 ; lAllard et all I2OOII : iHauschildt. Lowenthal. fc Baron 
2001f ). a generalized NLTE stellar atmosphere code, was 
utilized to compute the expected photospheric absorption 
line profiles for the primary of HD 44179. The calculated 
profiles included Nal. In § [5] we discuss the disagreement 
of the calculated line profiles with the observed profiles. 
The parameters adopted for this calculation were Te// = 
7,700 K, log g = 1.10 , [Z/H] = -3.5, [CNO/H] = -0.5 
fW aelkens et al. jl996l ). and a stellar mass of 0.8 Mq 



(Wit t et al.l 120091 ). The abundance determination utilized 
several Nal 8194.824 A line profiles with abundances from 
—0.10 dex to —1.00 dex. These lines were generated with 
the same atmospheric parameters, with the exception of 
their variable sodium abundance. Using the measured 
equivalent width of the Nal 8194.824 A line, a linear 
interpolation of the equivalent width values from the model 
lines was used to determine the Na abundance to be [Na/H] 
= -0.18 ±0.05 dex. 



Sodium is not significantly depleted in the spectra of 
HD 44179, while Fe and other refractory elements are sub- 
stantially depleted. The peculiar abundances of HD 44179 
were pointed out by Ivan Winckel. Waelkens. fc WatersI 
(|l995h . In our data there is no evidence of interstellar sodium 
in the line-of-sight to HD 44179. The observed depletion 
pattern i n HD 44179 i s similar to that of the interstellar 
medi u m (IMortonl 1974: Maa s. van Winckel. fc Lloyd EvansI 
I2OO5I : iGielen et al.1 :2009. : .Giridhar et al.l I2OI0I ). Our ob- 
servation of the Na abundance is consistent with the 
extremely metal- p oor p os t-AGB stars as di scusse d by 
Mathis fc LamersI (tl993). [Mathis fc Lamerd (|l992l ) and 



Waters. Trams, fc WaelkensI (| 19921 ) have proposed a mech- 



anism in which a circumbinary disc provides the most 
favourable conditions for 'cleaning' of the atmosphere of the 
refractory elements. In their model, refractory elements con- 
dense into dust and are blown away by radiation pressure. 
The now refractory-depleted gas is accreted back on to the 
star, thus altering the atmospheric abundances. 
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Table 2. Orbital parameters 



Symbol 


Parameter_^ 


Ml 


M2 


Van Winckel95^ 


Units 


P 


period 


319.2 ± 0.2 


318.8 ± 0.6 


298 


days 


Vo 


system velocity 


18.6 ± 0.1 


18.4 ± 0.5 


21.5 ± 1.2 


km s"-"^ 


Ki 


semi-amplitude 


9.8 ± 0.2 


14.9 ± 0.8 


13.0 ± 1.7 


km s"-*- 


e 


eccentricity 


0.276 ± 0.014 


0.320 ± 0.041 


0.45 ± 0.13 




Ul 


longitude of periastron^ 


-20.1 ± 3.2 


-f 7.3 ± 7.7 


28 ± 16 


o 


aism(i) 


semi-major axis 


0.281 ± 0.006 


0.414 ± 0.081 


0.32 


AU 


f(m) 


mass function 


0.028 ± 0.002 


0.049 ± 0.049 


0.049 


Mq 


ai 


true semi-major axis 


0.490 ± 0.006 


0.722 ± 0.081 




AU 


Mt 


total system mass 


1.49 ± 0.06 


1.80 ± 0.37 




M© 


m2 


secondary mass 


0.69 ± 0.04 


1.00 ± 0.37 




Mq 



'^ Orbital parameters for HD 44179 calculated with VCURVE, uncertainties are generated by the program except for the 
last five parameters. The uncertainties from VCURVE are used to compute the uncertainties in the last five parameters. The 
period of 317 days was used as a starting parameter for the fit to the phase folded data. Column Ml (method 1) and M2 
(method 2) are for different radial velocity determination methods discussed in detail in § [S] 
^ Ivan Winckel. Waelkens. fc Water3 l ll995^ 
'^ Referenced to periastron. 



5 PHOTOSPHERIC LINE SHAPES AND 
LIGHT PROCESSING 

5.1 Line measurements and asymmetries 

iHobbs et al.l ([2004) noted that the unsaturated photospheric 
absorption lines for all species have a similar FWHM of 
about 38 km s~^. After analysing our entire dataset we no- 
tice that not only are the FWHM for all absorption lines 
similar, the orbital phase dependence of the absorption line 
asymmetry is also the same; see Fig.[5]for examples. During 
the orbital phases of superior and inferior conjunction the 
photospheric line profiles are symmetric, while during the 
phases of maximum and minimum red-shift the line profiles 
are the most asymmetric. The Nal 8194.824 A line has a 
FWHM of 35.6 ± 4.1 km s"^ and the Ci 5052.167 A has 
a FWHM of 36.9 ± 1.3 km s~\ both in agreement with 
the previously measured values of iHobbs et al.l (|2004l ). The 
equivalent width of the Nal 8194.824 A line is 80.7 ± 3.0 mA 
and the equivalent width of the Cl 5052.167 A line is 
167.0 ± 5.7 mA. Neither shows significant dependence on 
orbital phase. While the FWHM and equivalent width show 
no orbital phase dependence, the line profile shapes (asym- 
metries) are phase dependent (Fig. [S]). In the following sub- 
section we will discuss various mechanisms, which may or 
may not shape the photospheric absorption lines. 



5.2 Broadening mechanisms 

The PHOENIX code produces lines that have the correct 
equivalent width, but are much narrower than the observed 
38 km s^^. Lines this broad are not unique among F-type 
post-AGB stars. The model lines do not include any broad- 
ening mechanisms other than thermal and pressure broad- 
ening. Significant rotation is not expected for a giant star. 
While it is possible the star is rotating, rotation cannot ex- 
plain the orbital phase-dependent asymmetry of the lines, 
as rotational broadening is symmetric. 

Macroturbulence is known to produce lines with sym- 
metric cores and asymmetric wings that conserve equivalent 



width; however, the required macroturbulent velocity would 
need to be roughly 23 km s~^, which would be very su- 
personic for this star's atmosphere. The isothermal sound 
speed is approximately 8 km s^^. Macroturbulence alone 
cannot explain the dependence of the asymmetry on the or- 
bital period. It also does not explain the extent of the wing 
either 40 km s~^ blue- or red-shifted depending on the or- 
bital phase. 

The primary star is only observable via the indirect line- 
of-sight. The location and nature of the scattering material 
is uncertain. If the material from which the light is scattered 
has some intrinsic velocity dispersion, then the observed line 
profile will be broadened as it scatters off this material. The 
scattering from the north and south of the nebula likely 
occur at slightly different scattering angles due to the tilt 
of the nebula with respect to the line-of-sight. However, the 
scattering from each end of the nebula likely occurs over a 
range of angles of order 5° . The fact that the binary stars are 
in motion around the centre-of-mass implies that the mean 
scattering angle with change slightly over the course of the 
orbit. In order to properly explore the scattering in the RR, 
a 3-D Monte Carlo model is planned to further study this 
system. 



5.3 Orbital phase-dependent asymmetries 

A double-lined spectroscopic binary could produce the line 
asymmetries we observe. The weaker component would be 
due to the companion star, while the stronger component 
would be due to the primary star. However, it is unlikely 
that the weaker component is due to absorption in the atmo- 
sphere of the companion star for the following r easons. The 
comp anion star is a 0.9 M© main sequence dwarf (|Witt et al.l 
|2009|). This implies that its luminosity is of order L©, while 
the primary star has a luminosity of order 6,000 L©. The 
ratio of the luminosities would mean an extremely weak con- 
tribution to the line from the secondary, while the observed 
contribution of the weaker component to the line is between 
1/3 and 1/2 of the flux of the strong component. It is also 
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unlikely that the weak component can be attributed to ab- 
sorption in the accretion disc around the secondary. While 
the accretion disc has a luminosity around 300 Lq , the con- 
tribution would still be smaller than 1/3 of the strong com- 
ponent. In addition, the m aximum effective temperature of 
the disc is about 23,000 K (|Witt et al.ll2009l '). and one would 
not expect atmospheric lines from the disc to have the same 
ionization balance as the atmosphere of the primary, which 
has an effective temperature of 7,700 K. Yet, observations of 
various photospheric absorption lines at different ionization 
states and species all exhibit the same asymmetric profiles. 

Pulsations in the atmospheres of R Coronae Borealis 
(RGB) stars show a second absorption-line component that 
is shifted toward either longer or shorter wavelengths of 
the centre of the photospheric absorption l ine, producing a 
variable asymmetric absorption-line profile (JRao fc Lambert! 
120081 ). These asymmetric profiles are reminiscent of the ob- 
served absorption-line profiles observed in HD 44179. How- 
ever, HD 44179 is not an RGB star. The brightness of RGB 
stars can decrease by nearly 8 mag at unpredictable in- 
tervals. HD 44179 has a very periodic and predictable be- 
haviour and does not experience such drastic changes in 
brightness. The brightness of HD 44179 changes by rou ghly 
0.14 mag (iWaelkens et al.ii 19961 ). iRao fc LambertI |200S ') at- 
tribute the line doubling observed in RGB stars to shocks 
in their atmospheres. Such atmospheric shocks are not lim- 
ited to RGB stars, similar line splitting is also seen in RV 
Tauri stars. Line splitting is also observed in other stars 
with loosely bound atmospheres. Gomplex atmospheric mo- 
tions are see n in other F-type post-AGB stars, for example 
SAO 96709 JLebre et al.lll996t ). 

We interpret the asymmetric wings as motion in the 
atmosphere of HD 44179. Material in the atmosphere is sys- 
tematically flowing toward the primary during some phases 
of the orbit and away from the primary during other phases. 
When the primary moves through the phase of periastron, 
the asymmetry is always on the blue side of the line core. 
The asymmetry is on the red side of the line core, when the 
primary is moving through the phase of apastron, compare 
Fig. [S] and Fig. [2] There is no discernible asymmetry at the 
conjunctions. The direction of flow, relative to the primary, 
appears to be affected by the varying size of the Roche lobe, 
see panel a of Fig. [TD] The Roche lobe radius of the primary 
at periastron is 0.34 AU. The binary components' separation 
at periastron is 0.91 AU. At apastron the stars are separated 
by 1.8 AU and the Roche radius of the primary is 0.66 AU. 
While the model atmosphere radius of the star does not 
exceed the Roche lobe radius, it should be noted that radi- 
ation pressure c an significantly redu c e the size of the effec - 
tive Roche lobe JDermine et ai.ll2009h . lDermine et al.1 (|2009l ') 
showed that the radiation pressure effect becomes stronger 
as a giant star's effective temperature increases. The authors 
only compute models for Te// < 5,000 K, but HD 44179 
has a higher Te//. Therefore, it is likely that the effective 
Roche lobe radius is decreased by this effect. The size of the 
Roche lobe can also be posi tively or negatively affected by 
the r otation rate of the star (jSepinskv. Willems. fc Kalogeral 
120071 ) in the case of systems with non-synchronously rotat- 
ing stars in eccentric orbits, as is the case with HD 44179. 
The rotation rate of the primary star is unknown; however, 
significant rotation is not expected in this system. 

Enhanced mass-loss occurs near periastron passage, but 



some of the lost material can be back-accreted onto the 
mass-losing star within that half the orbital period when 
the Roche lobe reaches its largest extent, increasing its vol- 
ume about eight-fold compared to its size at periastron. The 
self-accret ion is likely c aused by the ex panding Roche lobe 
(Scpinskv et al.ll201(]| ). ISepinskv et al.l (2010) showed that 
for systems with initial eccentricities greater than 0.2 and 
a mass ratio greater than 0.006, mass overfiow always leads 
to self-accretion. This suggests that we may be seeing the 
self-accretion occurring in HD 44179. In panel b of Fig. [10] 
we plot the derivative of the Roche lobe radius with respect 
to time. Negative slopes indicate a contracting Roche lobe, 
while a positive slope indicates an expanding Roche lobe. In 
panel c of Fig. [10] we show the velocity of the weak compo- 
nent of the asymmetric line, with respect to the star's sur- 
face. Negative velocities indicate material moving away from 
the star, while positive velocities indicate material moving 
toward the star. 



5.4 Light processing 

A final consideration is the uncertainty in the effective in- 
clination angle. While it is clear that the optically thick 
circumbinary disc prevents a direct view of the binary, the 
exact path or paths though which we observe the star affect 
our interpretation of the system. The cartoon model pre- 
sented in Fig. [l] and in other works previously cited, con- 
tains a rather thick, fiat circumbin ary disc. Such discs h ave 
been seen for example in HR 4049 (|Dominik et al.ll200l ). If 
the disc is flared, we could likely see into the system with a 
much higher effective inclination angle. The detailed geome- 
try of the inner region of the RR circumbinary disc requires 
a 3D model that includes resonant line scattering, which is 
beyond the scope of the present paper. 



5.5 Technique of photospheric line removal 

A specific model for how the RR produces these orbital 
phase-dependent photospheric lines is not needed in order 
to remove them from the spectra. Therefore, if we utilize a 
mostly empirical method to remove the absorption line com- 
ponent from the observed Nal D-line profile, we will be left 
with the nebular contributions to the spectrum. 

Using Fourier analysis, we can construct a convolution 
kernel that will take model atmosphere lines produced by 
PHOENIX and modify them to look like the observed absorp- 
tion lines in the spectra of HD 44179. Under the assumption 
that the model line profiles produced by phoenix have the 
correct shape for the underlying intrinsic photosphere, the 
PHOENIX model Nal D-line was taken as the 'true' shape 
of the line. A representative photospheric absorption line is 
then used as a template of the line shape for a given obser- 
vation. The observed G I 4932.049 A served as the template 
line shape. The observed G I 4932.049 A line was chosen 
because it is not contaminated by telluric features and it 
is unsaturated. We were unable to utilize the IR Nal line 
due to contamination of the line by telluric features during 
several orbital phases. 

For each observation, the phoenix model G I 4932.049 A 
line was Doppler-shifted to the appropriate wavelength ac- 
cording to the measured strong component of the RV. A 're- 
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Figure 11. Shown are the various steps in the model-subtraction 
process for the 5889.950 A line. This representative spectrum, 
# 15, was taken at maximum red-shift of the primary, which oc- 
curs at phase <p = 0.423, see Fig. [2] The normalized flux is in 
arbitrary units. The dotted line is the photospheric model line 
calculated by the PHOENIX code. The dashed line is the recon- 
structed photospheric line, which is the PHOENIX line multiplied 
by the response kernel constructed for the RR via Equation [T] 
The line with +'s is the observed spectrum, which consists of two 
averaged spectra with a total exposure time of 1 hour. The con- 
tinuum S/N is approximately 300. The solid line is the resulting 
model-subtracted profile. 



sponse' kernel was then constructed for the RR. The Fourier 
transform of the response was generated using Equation [T] 

-F(Oci) 



T{R) 



(1) 



-F(Mci) ' 

where R is the response, Oc i is the observed Ci profile, 
and Mc I is the model Cl profile. The Fourier transform of 
the Nal D-line model profiles were then multiplied by J- (R) 
resulting in an observationally generated reconstruction of 
the photospheric Nal D-line absorption profile. Once the 
reconstructed line is generated it can then be subtracted 
from the observed Nal D-line profiles, thus removing the 
effect of the underlying stellar absorption. An example of 
the PHOENIX model Nal D-line, reconstructed line profile, 
and subtraction are shown in Fig. [TT] 



6 SODIUM D-LINE PROFILES 

The Nal D-line profiles are a complex and variable superpo- 
sition of absorption and emission components, which repeat 
with orbital period. In Fig. [4] both observed Nal D-lines are 
shown at two different orbital phases to demonstrate the ex- 
tremes in the variation of the line profiles before the removal 
of the photospheric component. In Fig.[6l the high-resolution 
spectra taken at McDonald Observatory shows how distinct 
the various components of the line are. The complex ob- 
served profile is a superposition of: (1) a blue-shifted (with 
respect to the centre-of-mass RV) absorption from the mass- 
loss of the primary; (2) the photospheric absorption line; (3) 
an emission from the northern lobe of the bi-conical nebula; 
(4) and an emission from the southern lobe of the bi-conical 



nebula. The details of the various components and their or- 
bital variation will be discussed in the following subsections. 

6.1 Absorption components 

6.1.1 Red-shifted absorption 

Perhaps the most intriguing part of the Nal D-line profile 
(Fig. [4| is the red-shifted absorption that first appears in 
observations near inferior conjunction {(p = 0.21) and dis- 
appears near superior conjunction {(f> = 0.59). One possible 
explanation for this absorption feature is infalling material 
along the line-of-sight to the star that is present for only 
part of the orbit. If material is falling back onto the star 
as described in § 15.31 it has the wrong phase dependence to 
explain the red-shifted absorption. When the red-shifted ab- 
sorption is present, material is actually fiowing away from 
the star (Fig. llOf) . A more likely explanation, therefore, is 
the Doppler shifted photospheric absorption line. When the 
discrepancy between our phoenix model line FWHM and 
the observed FWHM is taken into account, as discussed in 
§ [5] the stellar absorption line seems the most likely expla- 
nation for the red-shifted absorption feature. After applying 
the subtraction technique in ij |5.5l the effect of the underly- 
ing photospheric absorption line can be removed, as shown 
in Fig. 111! The process completely removes the appearance 
of the red-shifted absorption component of the Nal D-line 
profiles for the entire range of orbital phase for which it ap- 
peared. As a result, the infalling material explanation can 
be ruled out, thus leaving the Doppler shifted stellar ab- 
sorption line as the explanation for the periodic red-shifted 
absorption feature. However, infalling material still remains 
the most likely explanation of the red-winged asymmetry 
seen in the photospheric absorption lines during the phases 
interval 0.58 < d < 1.2. 



6.1.2 Blue-shifted absorption 

The blue-shifted Nal D-line absorption is present at all or- 
bital phases, as shown in Fig. |4l The equivalent width does 
change as a function of orbital phase, but it is never zero, as 
seen in Fig. 1121 The equivalent width is smallest near phase 
(f> — 0.1, with an equivalent width of about 0.04 A. The ab- 
sorption is greatest near superior conjunction (</) = 0.59), 
with an equivalent width of approximately 0.085 A. The 
time during which the equivalent width is increasing strongly 
coincides with the observed periodic net outfiow from the 
primary star, as shown in panel c of Fig. 1101 The material is 
seen in absorption against the photosphere of the primary. 
Therefore, the absorbing gas must be in the indirect line- 
of-sight to the primary star, suggesting that the material is 
located in the cavity of the circumbinary-binary disc, not 
confined to the orbital plane. 

Another constraint on the physical location of the gas 
is the peak velocity of the blue-shifted absorption. Rela- 
tive to the system's centre-of-mass, the velocity of the blue 
absorption maximum is —20.7 ± 0.4 km s~^ as measured 
in the APO data and —21.4 ± 0.1 km s~^ as measured in 
the McDonald data. No dependence on the orbital phase is 
observed. The terminal velocity of the blue-shifted absorp- 
tion is near —50 km s"'^. This implies that the gas seen in 
the blue-shifted absorption against the primary star is not 
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Figure 12. Blue-shifted absorption equivalent width as a func- 
tion of phase measured in the model-subtracted 5889.950 A line. 
The uncertainty in the measurement is estimated to be about 
0.01 A. The vertical lines represent the phases of inferior conjunc- 
tion (IC) = 0.21, periastron (P) = 0.29, superior conjunction 
(SC) <f> = 0.59 and apastron (A) cj> = 0.79. Fig. [2] illustrates the 
key orbital locations. 



co-moving with the primary star, but rather moving at a 
fixed velocity away from the bary centre of the system, to- 
ward the observer along the line-of-sight. Therefore, the gas 
seen in absorption is outside the Roche lobe of the primary 
and is likely gas in the general outflow. We attribute the 
blue-shifted absorption to the general mass-loss driven out- 
flow from the primary star within the central cavity of the 
circumbinary disk. A lower limit of the total mass-loss rate 
through the bipolar outflow is given by the mass accretion 
rate of the disc around the secondary. The implied accretion 
rate (2-5xlO~^ M© yr~^), and therefore lower limit to the 
mass loss, is constrained by the required Lyman-continuum 
necessary to produce the obs erve d free-free r a dio em ission 
l|jura. Turner, fc Balml (|l997h and I Witt etHI (|2009l) ). 



6.2 Emission components 

6.2.1 Velocities 

Fig. U reveals double-peaked although not fully resolved 
emission components in the Nal D lines. The peaks are 
separated by about 12 km s~^ and are nearly symmetric 
about the system velocity. This is confirmed by the higher- 
resolution spectra shown in Fig. [B] In the centre-of-mass 
frame the velocities remain constant with orbital phase. This 
indicates that this gas is not involved in the binary motions 
of the star. The blue-shifted component is Doppler shifted by 
about —6 km s~^ with respect to the centre-of-mass frame, 
while the red-shifted component is Doppler shifted by about 
-f6 km s~^. These double-peaked emission lines were pre- 
vious ly attributed to m aterial in an edge-on circumbinary 
disc IjHobbs et al.ll2004l ). Material inside the cavity of the 
optically-thick circumbinary disc can only be viewed via 
the indirect line-of-sight scattering geometry proposed by 



IWaelkens et al.l (|l996l ). Viewing a disc at the effective incli- 
nation angle of 35° , one would expect to see a single broad- 
ened emission line from a disc, not a double-peaked profile. 
We should point out that the 12 km s~^ full width of the 
CO lines observed bv^ Jura. Balm, fc Kahand (|l995l ) is not 
due to the bipolar outfiow; the spatial distribution of the 
CO e mission confirms that i t originates in the circumbinary 
disc jBuiarrabal et al.ll2005l ). 

Close examination of the McDonald spectra. Fig. [6l re- 
veals that the narrow emission features sit on top of the pho- 
tospheric absorption line. The central absorption goes below 
the continuum, which is consistent with the location of the 
photospheric absorption line during that orbital phase {(p = 
0.814). The blue-shifted absorption is clearly separated from 
the photospheric line and the emission components. These 
data also show that the emission components have shoul- 
ders. The main emission peaks are separated by 11.6 km s^^. 
The blue-shifted emission peak occurs at —4.6 ± 0.2 km s^^ 
with respect to the centre-of-mass frame. The shoulder on 
the blue emission occurs at -7.3 ± 0.4 km s"\ The FWHM 
of the blue-shifted emission is 3.2 km s~^. The red-shifted 
emission occurs at +4.3 ± 0.2 km s~^ with the shoulder 
occurring at +7.7 ± 0.5 km s"\ The FWHM of the red- 
shifted emission is 2. 6 km s~^ and excludes the shoulder. 
IWhiteoak fc Gardneij ([1983) present a single averaged spec- 
trum of the two Nal D-lines at R = 120, 000; their results 
are consistent with our higher-resolution spectra. 

Since the emission line centres are nearly symmetric 
about the centre-of-mass RV, this suggests the possibility 
that we are seeing, simultaneously, outflowing material mov- 
ing both toward and away from the observer. The Nal D- 
lines are easily exited resonance lines that arise in material 
situated well past the inner cavity of the disc. The coro- 
n agraphic, long-slit Nal D STIS spectrum and the work 
of ISchmidt fc WittI l|l99ll ) observationally conflrm that the 
Nal D-line emission extends to large distances from the cen- 
tral source of the RR. From the STIS spectrum, the emis- 
sion can be traced to about 3.5" (at 710 pc this is roughly 
2,500 times larger than the binary orbit) n orth and south of 
the ce ntral source. The long-slit spectra of ISchmidt fc Witti 
(' 199l| ) reveal Nal D emission as far as 10" north and south 
of the central object. In addition, the STIS Nal D-line spec- 
trum shows blue-shifted emission emerging only from the 
southern part of the RR, while the red-shifted emission only 
originates from the northern part, as shown in Fig. [T] The 
northern view shows a weak absorption component on the 
blue side of the red-shifted emission close to the central 
source but not at larger distances from the optically-thick 
circumbinary disc. This is consistent with outflowing mate- 
rial seen against the star in the inner regions of the nebula, 
while at larger distances it is only seen in emission. 

Current ground-based observations of the RR lack the 
spatial resolution of STIS, thus incorporating light from 
both the north and south lobes of the nebula in a single 
composite spectral profile. The STIS Nal D-line spectrum 
demonstrates that the double-peaked emission profile we ob- 
serve from the ground originates from the extended biconi- 
cal nebula. The material in the bi-cones originates from the 
mass-loss driven outfiow of the primary star. The same ex- 
planation may also extend to all th e other double - peake d op- 
tical emission profiles reported bv iHobbs et al.l (|2004l ). but 
this has not been directly confirmed in the present work. 
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The double-peaked emission clearly does not originate from 
material orbiting in a c ircum binary disc as previously sug- 
gested bv lHobbs etHI \2QQ4 ). 



6.2.2 Inclination angle 

The range in inclination angles (89° to 79°: 90° can be 
safely ruled out in light of the double-peaked emission line) 
leads to a range in bipolar outflow speeds from about 300 
to 30 k m s~^. The best-f i tting value of the inclination an- 
gle from lBuiarrabal et al.1 (|2005l ). 85°, corresponds to about 
70 km s~^, which is of the order of the estimated escape ve- 
locity from the primary star, 83 km s~^. If we assume that 
the outflow has a velocity equal to the escape speed from the 
surface of the star, then we can infer the inclination angle 
to be 86°. Therefore th e value for the inclination used by 
iBuiarrabal et alj (|2005l ) would seem to fit our data best as 
well. The range of velocities from 30 km s"'^ to 70 km s~^, if 
interpreted as escape velocities, would correspond to a range 
in distances from the centre of the primary star of 1.6 AU 
to 0.3 AU respectively. This would seem to suggest that 
the blue-shifted absorption, with a maximum velocity near 
50 km s~^ fij I6.1.2|) . is produced by the same outflow that 
produces the emission lines. Therefore, with the Nal D-lines 
we observe the mass-loss driven outflow, powered by the pri- 
mary star, from inside the cavity of the optically-thick cir- 
cumbinary disc out as far as 10". Inside the cavity we see 
the material in absorption against the primary star, while 
outside the cavity we directly view the material in emission. 



6.2.3 Emission strength 

The height of the centroids of the blue-shifted and red- 
shifted emission peaks above the continuum for both lines 
in the Na I D doublet were measured as a function of phase 
for our model-subtracted spectra. The height of the blue- 
emission peak is referred to as the V-height, while the red- 
shifted height is called the R-height (both are plotted in 
Fig.[T3]as a function of orbital phase). These quantities allow 
us to examine the components of the emission proflle sepa- 
rately and they are preferable to using the equivalent width. 
The equivalent width could in principle be used to quantify 
these components. However, since the emission peaks are 
blended and there are absorption features on both sides of 
the emission peaks it complicates the interpretation of the 
equivalent width measurements. The equivalent width for 
the blue-shifted and red-shifted emission components can 
only be determined through deblending, which increases the 
uncertainty in the measurements. This is due, primarily, to 
the uncertainty of the underlying shape of each emission 
component. 

The outflow, as characterized by the V-height (FiE. ll3|) . 
is strongest near the phase of maximum RV {(j> = 0.4). The 
R-height is greatest near phase = 0.6. The ratio of the 
V-height to the R-height (V/R-ratio) was also used to char- 
acterize the variability. Fig. [14] shows the V/R-ratio for the 
5890 A D2-line, which peaks near phase </> = 0.4. The V/R- 
ratio for both the STIS and McDonald data are included 
in Fig. [14] and show good agreement with the data from 
APO. In order to make a comparison with the APO data, 
the McDonald datum was photosphere-subtracted and arti- 
ficially degraded to match the resolution of the APO data. 
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Figure 13. Plotted are the V-height (squares) and the R-hcight 
(triangles) for the 5889.950 A line measured in arbitrary units as 
a function of phase. The measurements have been corrected for 
the periodically-variable continuum level. These measured were 
performed on the model-subtracted spectra. The measurement 
error is on order of the size of the symbols. The scatter in the 
data may represent small time scale variations and uncertainty 
in the continuum normalization of the data. The vertical lines 
have the same meaning as in Fig. 1121 Fig. [2] illustrates the orbital 
locations. 



The STIS data are spatially resolved and the contribution 
from the stellar atmosphere should be negligible due to the 
choice of extraction apertures. The variability of the V/R- 
ratio is dominated at all phases by the blue-shifted compo- 
nent. Therefore, it is not surprising that the overall phase 
dependence is observed to be the same, peaking near phase 
(j) = 0.4. 

The ratio of the flux between the D2-hne (5890 A) and 
the Dl-line (5896 A) for optically thin processes is the ra- 
tio of the oscillator strengths, which is equal to two. For 
optically-thick processes the ratio approaches unity. The ra- 
tio of the V-height 5890 A to V-height 5896 A is measured to 
be unity and shows no phase dependence. This implies that 
the emission must be optically thick. The blue absorption 
equivalent width (W), however, does obey the optically thin 
limit of W5890 / W5896 = 2, which would seem to disagree 
with the emission ratio. In the McDonald data the W ratio is 
1.4 at its higher-resolution. When the McDonald data were 
artificially degraded for comparison, the W ratio increased 
to 1.5. This suggests that in the lower resolution APO data 
the blue-shifted emission is affecting the W ratio and that 
the material seen against the star is not optically thin, but 
not entirely optically thick. Taking the McDonald data into 
consideration, it seems that we are seeing the same outflow 
of material in the cav i ty tha t is flowing into the outer neb- 
ula. ISchmidt fc Witt! (|l99ll 'l observed that the Nal D-line 
emission was strongest along the bi-cone walls or 'whiskers'. 
At large distances from the centre of the nebula we directly 
see the emission, which may be dominated by emission from 
walls of the outflow cavity. The optical depth through the 
edge of outflow would be largest. These observations are con- 
sisten t with the work of iWarren-Smith. Scarrott. fc MurdinI 
(jl98ll ) who found that the sodium emission was optically 
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thin on either side of the 'whiskers' and opticaUy-thick on 
the 'whiskers'. 

Similarly, the ratio of red-shifted emission height R- 
height 5890 A to R-height 5896 A is also constant with 
phase; however, the ratio is about 0.84. The 5896 A Dl-line 
consistently has a stronger red component than the 5890 A 
D2-line. This remains a puzzle; it was expected that the 
red-shifted emission component may be more highly scat- 
tered and thus weaker than the blue-shifted component. 
This effect should apply equally to both lines of the doublet. 
Since the ratio for the blue-shifted emission is unity, imply- 
ing near optically-thick emission, the expectation would be 
that the red emission ratio would also be unity. However, 
the smaller ratio of the red component is unexpected for 
the optically thick case. This could indicate that the pho- 
tospheric subtraction technique of § 15.51 is insufficient, but 
without the means to spatially separate the components at 
high-resolution we cannot resolve the issue in this paper. 

The blue-shifted emission and blue-shifted absorption 
in the APO spectra are seen to be strongest near phase <j) = 
0.48 (just after periastron) coinciding with the period of the 
fastest mass outflow from the primary star, based on the 
asymmetry of the photospheric absorption lines, see panel 
c of Fig. [To] During these phases, mass is being added to 
the outflow resulting in increased absorption and emission. 
The close proximity in time of the strongest emission and 
absorption is suggestive that the outflows seen in absorption 
and emission are one and the same. This argument provides 
further support to the discussion in ij 16.21 relating the veloc- 
ities of these components. Note that the gas seen in Nal D 
absorption is signiflcantly slower than that of Ha, seen in 
absorption near superi or conjunction, a nd attributed to the 
jet from the secondary (|Witt et al.ll2009l V This also supports 
the idea that the outflow seen in Na I D emission and absorp- 
tion originates from the primary star, rather than from the 
secondary star's accretion disc jet. However, why the emis- 
sion components, which extend to large distances from the 
central source, are correlated with the motion of the primary 
star is unknown. 



7 SINGLE VS. DOUBLE PEAKED EMISSION 

A comparison of single-peaked emission (for example Ha 
IjWitt et al.ll2009l ')) and double-peaked emission (as seen in 
Nal D-line) reveals that they are representative of differ- 
ent environments in the RR system. The STIS Ha emis- 
sion, seen in Fig. [S] can be traced to roughly 1.3" north 
and south of the central source. The Ha has a smaller spa- 
tial extent and is seen primarily through scattering (indi- 
rect line-of-sight). The Nal D-line emission can be traced to 
3.5" north and south of the central source in the STIS spec- 
tra. The large spatial extent of the Nal D-lines allows them 
to be seen both directly (at larger distances from the cen- 
tral source) and indirectly via scattering. The single-peaked 
Ha emission is only seen indirectly and, therefore, is arising 
only from the i nner regions of the RR. The single-peaked 
emission fines of lHobbs etakl (|2004l ) and lWitt et akl (|2009l ') 
most likely originate f r om th e Hll region as mentioned by 
IJura. Turner, fc Balml lll997f). 

In table 2 of lHobbs et al.l (|2004l 'l they classify 84 emis- 
sion fines from 12 atomic species based on the shape of their 
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Figure 14. V/R ratio for the 5889.950 A line in arbitrary units 
as a function of phase, measured from the model-subtracted spec- 
tra. The scatter seen in this data as shown in Fig. 1131 disappears 
because this is a relative measurement on each spectrum. Small 
fluctuations on short time scales, and uncertainty in the nor- 
malization of the spectra are removed. The measurement errors 
are of the order of the symbol size. The vertical lines have the 
same meaning as in Fig. 1121 The solid squares are from the APO 
dataset. The open triangle is from the artificially broadened Mc- 
Donald spectrum. The open inverted triangle is from STIS spec- 
trum. Fig. [2] illustrates the orbital locations. 



profile. They list 24 'narrow', 12 'broad' and 48 'double- 
peaked' atomic lines. The 'narrow' or single-peaked emis- 
sion spectra include the [Can] 7291.47 A lines. The Call 
3933.66 A K-line is also among the lines classified as single- 
peaked, though upon close inspectio n an asymmetry can 
be see n in the line. The spectra of IWhiteoak fc Gardneii 
|l983) reveal that, at higher- resolution, the Call K-line is 
in fact double-peaked with a separation between peaks sim- 
ilar to that of the Nal D-lines. Close examination of the 
double-peaked profiles reveals that the relative height of 
the blue-shifted and red-shifted emission peaks vary from 
species to species, for example the Nal D-lines and Can 
K-line. Most of t hese, however have near-equal fiux for the 
two com ponents |Hobbs et al . 2004). The Call K-line spec- 



tra of iWhiteoak fc Gardner! (1983. 1 shows a weaker blue- 
shifted emission component relative to the stronger red- 
shifted emission component, opposite to what is seen in Nal. 
The exact reason for this is unclear. The present work did 
not look at the orbital phase-dependence of the other emis- 
sion lines in the system. Future work will look at how the 
V-height and R-height ratios of other multiplet lines behave 
as a function of the orbital period. This could also be done 
with single lines in principle; however, with multiplet lines 
one assumes that each component line should be equally af- 
fected by scattering. One also knows what the fiux ratios 
between line components should be and also that the lines 
arise from the same atom; therefore, the environmental con- 
ditions giving rise to the lines are the same. Hence, any 
discrepancy in such measurements may reveal more about 
how this complex system behaves. 
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8 CONCLUSIONS 

The complex profiles of the Nal D-hnes in HD 44179, pe- 
riodically varying in intensity and shape with the orbital 
period as the photometric primary post-AGB component of 
the binary moves through its eccentric orbit, reveal detailed 
information about the varying mass-loss from the primary. 
Peak mass-loss occurs during the orbital phases around pe- 
riastron, to be replaced by partial re-accretion during the 
phases near apastron. 

The primary conclusions are: 

(i) The double-peaked sodium emission lines are pro- 
duced in the bipolar outflow, seen via the direct line-of- 
sight (inclination of 85°). The observed velocity difference 
of 12 km s~^ is a result of the 5° tilt of the bipolar outflow 
axis with respect to the plane of the sky. 

(ii) The same outflow seen in the double-peaked emission 
is also seen in absorption against the photospheric contin- 
uum of HD 44179 via the indirect line-of-sight (effective in- 
clination angle of 35°). The velocity of the outflow approach- 
ing the escape velocity from the primary suggests that the 
primary star is the ultimate source of the outflowing mate- 
rial. 

(iii) The single peaked emission lines, for example Ha, 
likely originate from the Hll region located in the cavity 
of the circumbinary disc. The material with single peaked 
emission lines is at rest with respect to the centre-of-mass. 

(iv) The red-shifted absorption component of the sodium 
profile, present for only part of the orbit near periastron, is 
most likely produced by the underlying photospheric Na I D 
absorption-line of HD 44179. The specific orbital behaviour 
is governed by the RV of the primary star. 

(v) The photospheric asymmetric blue-winged profiles in 
HD 44179 are consistent with enhanced mass-loss near pe- 
riastron passage. Some of the material lost is accreted onto 
the disc surrounding the secondary star. 

(vi) The photospheric asymmetric red- winged profiles are 
consistent with self-accretion of a portion of the material lost 
near periastron passage. This self-accretion occurs on a time 
scale of half the orbital period, centred around the phase of 
apastron. 

It is clear that the cause of the broadened photospheric 
absorption lines is of great importance for understanding 
and interpreting the spectra of HD 44179. With current ob- 
servations, some ambiguities cannot be resolved. The broad- 
ening mechanisms investigated in this paper are either un- 
likely or cannot fully explain the observed line-profiles. It 
is likely that there may be a superposition of mechanisms 
at play in shaping these spectral features. The indirect line- 
of-sight adds to the complication of interpreting the spec- 
tra. The effective inclination angle used here may in fact be 
somewhat different. It is more likely that we see HD 44179 
under a small range of angles, of which the effective inclina- 
tion angle is the mean. It is likely the case that the spread 
of viewing angles leads to broadened profiles. In addition, 
there is a slightly different viewing geometry for the north 
and south ends of the nebula. Spatially-resolved, higher- 
resolution spectra of the photospheric lines could help re- 
solve some of these remaining puzzles. A full 3-D model is 
called for in order to understand our complex view of the 
binary deep inside the circumbinary disc. 
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